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(57) Abstract: The invenrion concerns a method for real time navigation for locating a nomad (SUR) using radio signals with three 
carriers of three different frequencies, to determine the position of a user, transmitted by satellites (SAT r GPSE EJ A SAT, r GPSE EH ). 
The method comprises a first step which consists in determining extra- long-path carrier phase ambiguity, a second step which consists 
in estimating long-path phase ambiguity, and a third step which consists in resolving the phase ambiguity of one of the frequencies. 
One additional step consists in real time application of ionospheric corrections during the third step, said ionospheric corrections 
being based on a constantly updated ionospheric model of said ionospheric layer calculated by a reference fixed earth station (REF- 
REF E \ combined with geodesic data calculated by a so-called master reference fixed earth station {REFm-REF me ). The invention 
also concerns a system for implementing said method. 
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PROCED^ ET SYSTEME DE NAVIGATION EN TEMPS REEL A L'AIDE DE 
SIGNAUX RADIOELECTRIQUES A TROIS PORTEUSES EMIS PAR DES 
SATELLITES ET DE CORRECTIONS lONOSPH^RIQUES 

Uinvention concerne un procede et un systeme de navigation en 
temps reel a I'aide de signaux radioelectriques a trois porteuses 6mis par des 
satellites et de corrections ionospheriques, plus precis6ment de corrections 
obtenues a I'aide d'un modele ionospherique en temps reel, continument mis £ 
5 jour, le modele etant base sur des donnees d'un systeme de navigation par 
satellite, par exemple implements comme un modele de voxels en trois 
dimensions. 

Elle trouve une application particulfere, bien que non exhaustive, dans 
le domaine de la navigation instantanee tr6s precise, typiquement avec une 
10 precision inferieure au decimetre comme il le sera montr6 ci-apr£s, ce sur des 
distances de I'ordre de la centaines de kilometres ou plus. 

Pour fixer les idees, on se placera dans ce qui suit dans ce cas de 
Implication preferee de Pinvention, sans en limiter en quoi que ce soit sa 
portee. 

15 Une des techniques actuellement couramment utilisees pour obtenir 

un positionnement precis d'un "objet", fixe ou mobile, et, dans le dernier cas, 
des donnees relatives a sa position, son deplacement, la direction de ce 
deplacement, et/ou sa Vitesse, est d'utiliser des signaux radioelectriques ernis 
par des satellites artificiels toumant autour de la terre. Le terme "objet" doit etre 

20 compris dans son sens le plus general, notamment un vehicule terrestre, 
maritime ou a6rien. Pour simplifier, on appellera "nomade" cet "objet" dans ce 
qui suit. 

Diverses techniques sont connues pour obtenir le positionnement 
precite. Elles sont notamment basSes sur la connaissance de la position 
25 instantan6e de plusieurs satellites dans Tespace (ou constellation de satellites, 
ces satellites pouvant etre gSostationnaires ou non) et de la vitesse de 
propagation des ondes radioelectriques. Des horloges trSs precises sont 
embarqu6es a bord des satellites et les signaux emis comprennent des 
informations horodatees, ce qui permet de connaitre pr6cis6ment Sgalement 
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I'instant demission et de reception. II est alors possible de determiner la 
distance th6orique separant un satellite en vue du nomade d un instant donn6 
de celui-ci, connaissant la Vitesse de propagation des ondes et le temps mis 
pour atteindre le nomade. Si des satellites en nombre suffisant sont 

5 obseryables, on peut determiner ainsi les cordpnnees du nomade par rapport a 
un referentiel, en deux dimensions (longitude et latitude sur le globe terrestre), 
voire en trois dimensions (longitude, latitude et altitude : verticale). 

Cependant, comme il va I'etre montr6, du fait d'accumulation d'erreurs 
dans les mesures dues & diverses causes, les distances calcul^es ne sont 

10 qu'approximatives et la determination de la position du nomade souffre d'une 
imprecision plus ou moins importante, selon les technologies mises en oeuvre. 

Un des systemes de navigation par satellite les plus connu est le 
sysfeme dit "GPS" pour "Global Positioning system" ou "Systeme de 
Positionnement Mondial". 

15 Classiquement, les satellites "GPS" (ou d'un type semblable) 6mettent 

sur deux bandes de frequences, gen6ralement design6s Li (frequence de 
porteuse egale a 1,575 GHz) et L 2 (frequence de porteuse 6gale a 1 ,227 GHz), 
d'ou la qualification de "double frequence" ou "dual -frequency" qui leur est 
appiiqu&e. 

20 L'utilisation de ces deux frequences, conformement d certains 

precedes bien connus de I'Homme de M6tier, permet d'ameliorer la precision 
de la determination de la position d'un nomade par rapport a un referentiel, 
mais necessite des rScepteurs "GPS" plus complexes et plus onereux. 

Une determination de position peut etre faite, selon deux methodes 
25 principals : en temps r§el ou en temps differe, en effectuant ce qui est appel6 
un "post-traitement". Dans le premier cas, on parte egaiement de resolution 
"single-epoch" selon la terminologie anglo-saxonne la plus couramment utilisee 
(terme qui sera utilise ci-apr6s) ou de resolution instantanee, les calculs etant 
effectues lors d'une seule "6poque" tfobservation. La seconde methode ("post- 
50 traitement") permet d'ameliorer la precision. Cependant, si pour des nomades 
se deplagant lentement (un navire par exemple), cette derniere methode ne 
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presente pas d'inconvenients majeurs, elle n'est pas appropriee pour des 
nomades se d6plagant tres rapidement (aeronefs par exemple). 

On peut encore ameliorer la precision, en combinant les signaux emis 
par les satellites avec des signaux en provenance de stations terrestres fixes 

5 de reference, dont les positions sont parfaitement connues. Cependant, il est 
necessaire, si le nomade parcourt de grandes distances, que le reseau de ces 
stations soit relativement dense, toujours dans le cas ou une grande precision 
dans la determination de la position du nomade est souhaitable, ce qui 
augmente correlativement le coflt du systeme global. 

10 En outre, parmi les causes importantes d'erreurs, la refraction 

differentiate ionospherique, si on considere des distances 6gales & des 
dizaines de kilometres ou plus, est un des problemes principaux affectant les 
capacites de resolution instantanee d'ambiguTte de phase de porteuse, et par 
consequent la possibility de fournir une navigation dont la precision est de 

15 I'ordre du centimetre avec les systemes satellites de navigation mondiale a 
double frequence tels que le systeme "GPS" precite. Cette caracteristique 
restera vraie en ce qui conceme les futurs systemes mettant en oeuvre trois 
frequences, comme le systeme "GALILEO" et le systeme "GPS" modernise. 

En effet, les systemes de navigation satellites § trois porteuses en 

20 cours de planification offrent les avantages potentiels d'un taux eleve de 
reussite de la resolution instantan6e d'ambigurte et d'integrite dans celle-ci, 
avec une quantite de calculs geodesiques minimal e. Ceci est dG notamment au 
fait que Ton dispose simultanement d'un nombre plus important de donnees 
distinctes (tfest-a-dire associees aux trois frequences precitees), ce qui 

25 ameiiore correlativement les chances d'obtenir une resolution d'ambiguite 
instantanee ("single-epoch"). 

Mais, de nouveau, cette resolution peut etre serieusement affect6e par 
la refraction ionospherique, comme detailie ci-apres. 

Pour atteindre des precisions elevees dans la determination 

30 instantanee de la position d'un nomade, notamment d'un nomade se depla?ant 
sur de grandes distances, le besoin continue done £ se faire sentir de mettre 
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en oeuvre des techniques permettant, notamment, de reduire Finfluence 
nefaste de la refraction ionospherique. 

Dans I'etat de la technique, des precedes visant a repondre a ce 
besoin ont ete proposes. 

5 On peut citer notamment le proc6d6, dit "TCAR", de "resolution 

d'ambiguTte a trois porteuses". Ce proced6 est d6crit dans I'article intitute : 
"ANALYSIS OF THREE-CARRIER AMBIGUITY RESOLUTION (TCAR) 
TECHNIQUE FOR PRECISE RELATIVE POSITIONING IN GNSS-2", de U. 
VOLLATH et al., paru dans "Proceedings of the ION-GPS", 1998, IX-O-13, 

10 pages 1-6. 

On peut egalement citer le precede, dit "CIR", de "resolution de 
nombres entiers en cascade". Ce precede est decrit dans Particle de Jaewo 
JUNG et al., intitute "Optimization of Cascade Integer Resolutionwith three Civil 
GPS Frequencies", paru dans "Procedings of the ION-GPS" 2000. 

15 Ces deux techniques partagent une approche de base semblable : les 

ambiguites de double difference de nombres entiers sont successivement 
resolues en calculant des battements de frequences d'ondes. Ce calcul est 
effectue en partant du battement de la longueur d'onde la plus elevee a la plus 
courte, incluant des combinaisons de phases de porteuse de chemins dits 

20 "long" et "extra long" (avec des longueurs d'onde de 7,480 m et de 0,862 m, 
respectivement), et une premiere porteuse a la "frequence Li" (avec une 
longueur d'onde de 0,190 m). 

En particulier, le precede "TCAR" constitue une approche simple qui 
essaie de resoudre instantanement ("single-epoch") Tensemble complet des 

25 ambiguites. Mais les performances de "TCAR" sont fortement affectees par le 
decorrdlation ionospherique de refraction qui se produit avec la distance. En 
effet, comme details ci-apres, le retard d'ionosph§re est un probleme quand 
(comme d'ailleurs dans le cas des systemes a deux frequences) la valeur de sa 
double differentiate est superieure a 0,26 TECU (ce qui correspond a 4 cm de 

30 retard pour Li). 

"TECU" est une unit6 utilisee pour decrire certaines caracteristiques 
6lectriques de Tionosphere. En effet, I'ionosphere peut etre decrite a i'aide 
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d'une carte reprSsentant un compte de nombre total d'electrons, ou "TEC" 
(pour "Total Electron Content"). La carte represente Integration du nombre 
d'6lectrons suivant une direction verticale comme une fonction de la latitude et 
de la longitude. L'unite de 'TEC" est appetee 'TECU" (pour TEC's Units"), 

5 avec 1 TECU = 10 16 electrons contenus dans un cylindre aligns sur la ligne 
d'observation d'un satellite observe, dont la section est de 1 m 2 . Les particules 
chargees de I'ionosphere sont generees par le soleil, dont I'intensite du 
rayonnement varie naturellement en fonction de Pheure consid^ree. Du fait que 
la terre tourne autour de son axe, sous la couche ionospherique, la carte "TEC" 

10 est habituellement considSree comme representant une reference fixe par 
rapport au soleil, mais changeante en fonction de I'heure. 

Le seuil rappele ci-dessus peut etre facilement d§pass6 comme il est 
possible de le constater en consultant des cartes ionosph6riques ('TEC") de 
retards verticaux calcules a partir des donnSes "GPS". De telles cartes sont 

15 6ditees, a titre d'exemple par le "Jet Propulsion Laboratory", I'universite de 
Berne, etc., et publiees sur Internet par P "University Corporation for 
Atmospheric Research" et d'autres organismes semblables. 

En consequence, en vue d'amSliorer encore le procede "TCAR", une 
approche integree dite "ITCAR" (pour "Integrated TCAR") a ete developpee. 

20 Cette technique est decrite, par exemple, dans I'article de VOLLAT et al. 
pr6cite. 

Cette technique recourt a des algorithmes de recherche et a un filtre 
de navigation dans lesquels les ambiguites font partie des signaux de sortie et 
les erreurs r6siduelles ionospheriques sont grossierement estimSes. Pour une 

25 description plus detaillee des techniques mises en oeuvre, on se reportera avec 
profit a cet article. 

Cependant, bien qu'apportant une amelioration significative, la 
technique "ITCAR" reste n6anmoins affectee par le manque de connaissance 
de la double difference de refraction ionospherique, limitant ainsi le taux de 

30 reussite de la resolution d'ambiguYte pour des distances superieures a 
quelques dizaines de kilomdtre, comme d6crit dans I'article de VOLLATH et a I., 
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intitule : " Ambiguity Resolution using Three Carriers - Performance Analysing 
using "Real" Data", paru dans "GNSS Symposium", Seville, mai 2001. 

II a egalement ete propose, toujours dans le but d'am6liorer la 
precision de la determination de la position d'un nomade par rapport d un 
5 referentiel, de combiner un modele ionospherique en temps reel de 
I'ionosphere, obtenu a partir de donnees de deux frequences ("dual- 
frequency") gen6r6es par un reseau de stations fixes, avec des donnees d f un 
programme geodesique, et d'utiliser de telles donn6es pour effectuer des 
corrections ionospheriques. Ce procede a 6t6 utilise avec un certain succes 

10 pour resoudre des ambiguTtes en temps r6el dans les systemes a deux 
frequences du type "GPS". 

Un tel precede de ce type, dit "WARTK" (pour "Wide Area Real-Time 
Kinematics" ou cinematique en temps reel sur de larges etendues), est par 
exemple d6crit dans I'article Hernandez-Pajares et aL, intitule "Tomographic 

15 modeling of GNSS ionospheric corrections : Assesment and real-time 
applications", paru dans "ION GPS", 19-22 septembre 2000, pages 616-625. 
Ce precede permet, certes, d'attenuer les effets nefastes des perturbations 
subies par les ondes radioelectriques se propageant dans I'ionosphere et, par 
consequent, d'ameliorer tres sensiblement le taux de reussite de la resolution 

20 d'ambiguTte de phase et la precision de la determination de la position d'un 
nomade par rapport £ un referentiel, mais il necessite une quantite de calculs 
importante pour 6tabiir le module precit6 en temps reel et determiner, 
egalement en temps r6el, les corrections ionospheriques a apporter aux 
mesures de distances. 

25 [.'invention vise a pallier les inconvenients des procedes et systemes 

de Tart connu, et dont certains viennent d'etre rappetes. 

L'invention se fixe pour but un procede de navigation en temps r6el, 
tr6s precis (determination de la position d'un nomade typiquement meilleure 
que le d6cim6tre) d I'aide de signaux radioelectriques a trois porteuses emis 

30 par des satellites et de corrections ionospheriques obtenues a partir de 
donnees d'un modele de description ionospherique d'une region traversee par 
le rayonnement radioelectrique, combine a des donnees geod6siques. 
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Pour ce faire, selon une premiere caracteristique importante de 
l'invention, un precede simple de resolution instantanSe ("single-epoch") 
d'ambigurtes de phases de trois porteuses emise par des satellites, du type 
"TCAR" precite est utilise. Des corrections ionospheriques precises en temps 

s reel, basees sur un modele de description ionospherique, sont generees par 
des stations fixes de reference. Ces corrections ionospheriques sont 
transmises aux utilisateurs qui les incorporent a "TCAR". 
L'invention presente des avantages importants. 
Les dispositions, propres a l'invention, permettent d'obtenir de 

10 meilleures performances que le precede "ITCAR" sur de longues distances, £ 
une position tres 6loignee de la station de reference le plus proche. Le precede 
selon Tinvention permet typiquement la resolution instantanee de plus de 90 % 
des ambiguftes £ plus de 100 kilometres de cette station et pratiquement de 
100 % pour un eloignement de I'ordre de 60 km, ce dans un scenario de 

15 fonctionnement peu favorable : valeurs ionospheriques basses et conditions 
d'ensoleillement maximum. Les precedes de Tart connu mettant en ceuvre trois 
frequences permettent seulement, dans des conditions comparables, un taux 
de succes typiquement de I'ordre de 60 %. 

On congoit aisement que cette demfere caracteristique soit tres 

20 interessante, car elle ne necessite pas la mise en place d'un reseau de stations 
terrestres fixes de reference tres dense (stations eloignees les unes des autres 
de I'ordre de la centaine de kilometres). II s'ensuit un gain global tres important 
: dans le cas de I'Europe, par exemple, pour couvrir ce continent, il est possible 
d'economiser 99 % de stations de reference par rapport a un r6seau dont les 

25 stations de reference seraient seulement eloignees d'une distance de Tordre 
d'une dizaine de kilometres les unes des autres, tout en conservant une meme 
precision dans la determination de la position d'un nomade que celle obtenue 
par les precedes de Tart connu. C'est en effet cette densite d'implantation qui 
serait n6cessaire pour obtenir une precision de navigation meilleure que le 

30 decimetre, vis6e par le precede de l'invention, dans son application preferee. 
Enfin, le reseau de stations terrestres fixes de reference necessaire au 
precede de l'invention est tres similaire au reseau en cours d'implantation pour 
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le systeme europeen de navigation par satellites connu sous Pabreviation 
"EGNOS GNSS" (pour "European Geostationary Navigation Overlay System - 
Global Navigation Satellite Systems"). Ce systeme est base sur une 
constellation de satellites "IMMERSAT III", ce qui laisse penser qu'une 
5 compatibility importante existe avec les installations deja existantes ou en 
coursd'implantation. 

Par ailleurs, le proc6de selon invention reste tout & fait simple, 
comrne le procede "TCAR". En particulier, il n'exige qu'une quantite de calculs 
faible pour Putilisateur (recepteur du nomade), comparee a "ITCAR" et 
10 egalement a "WARTK", procedes developpes pour les systemes de 
positionnement a deux frequences, comme il a ete rappele ci-dessus. 

D'autre part, le modele de description ionospherique en temps reel de 
la region traversee par le rayonnement radioelectrique 6mis par les satellites 
observes utilise seulement des donnees de phase de porteuse, ces donnees 
15 etant combinees avec des estimations geodesiques calculees dans un centre 
unique, qui peut etre Pune des stations, drte "maitre", du r§seau de stations de 
reference. 

^invention a done pour objet principal un Procede de navigation en 
temps r§el § Paide de signaux radioelectriques a trois porteuses de premiere, 

20 deuxieme et troisieme frequences distinctes, de valeurs croissantes en allant 
de ladite premiere a ladite troisieme frequence, pour determiner la position d'un 
utilisateur, dit nomade, lesdits signaux radioelectriques etant emis par un 
nombre determine d'emetteurs embarques dans des satellites tournant autour 
du globe terrestre et en vue dudit nomade, lesdits signaux etant re$us par un 

25 recepteur associe au dit nomade et un r6cepteur associe a au moins une 
station terrestre parmi une plurality de stations terrestres fixes, dites de 
reference, lesdits signaux radioelectriques traversant une couche de 
Patmosphere dite ionospherique entourant ledit globe terrestre et y subissant 
des perturbations engendrant des ambiguTtes de phases desdites porteuses, 

30 caracterise en ce qu'il comprend au moins les etapes suivantes : 

une premiere etape consistant en la determination, dans ledit 
nomade, de PambiguTte de difference de phases entre lesdites troisieme et 
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deuxieme porteuses, drte de "chemin extra-long", a partir d'une combinaison de 
pseudo-distances utilisant une vaieur de code unique ; 

une deuxfeme etape consistant en I'estimation, dans ledit 
nomade, de I'ambiguTte de difference de phases entre lesdites premiere et 
5 deuxieme porteuses, dite de "chemin long", a partir de ladite ambiguite de 
"chemin extra-long" determinee lors de ladite premiere etape ; 

une troisieme etape consistant, dans ledit nomade, en la 
resolution de I'ambiguite d'une desdites frequence, a partir de ladite ambiguTte 
de "chemin long" estimee lors de ladite deuxieme etape ; et 
10 une etape supplemental I'application de corrections ionospheriques 

en temps reel pendant ladite troisieme etape, lesdites corrections 
ionosph6riques 6tant basees sur un module ionospherique en temps reel de 
ladite couche continument mis a jour. 

L'invention a encore pour objet un systeme de navigation par satellite 
15 mettant en ceuvre ce procede. 

[.Invention va maintenant etre decrite de fagon plus detaillee en se 
referant aux dessins annexes, parmi lesquels : 

la figure 1 A illustre schematiquement une coupe partielle du globe 
terrestre, de la couche ionospherique I'entourant, decoupee en voxels, 
20 ainsi que des faisceaux radioelectriques g6n6r6s par des emetteurs 

"GPS" de trois satellites de communication et trois stations terrestres 
recevant ces signaux ; 

la figure 1B illustre, de fa?on sch§matique, Tarchitecture d'un 
systeme de navigation complet, pour la mise en ceuvre du procede de 
25 l'invention, vue du c6te d'un recepteur "GPS" d'un nomade ; 

la figure 2 represente schematiquement la ligne de base et les 
positions d'une station terrestre de reference et d'un nomade, par 
rapport a des axes de coordonnees ; 

la figure 3 est un agrandissement (zoom) de la figure 2 montrant 
30 la trajectoire du nomade, par rapport a ces memes axes de 

coordonnees ; 
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la figure 4 est un graphique illustrant des exemples de retards 
estimes d'un modele ionospherique en temps reel selon I'invention 
compares £ des retards ionospheriques reels ; 

la figure 5 est un graphique illustrant des exemples de double 
5 difference de "Contenu Total d'Electrons Obliques" ("STEC"), estimee 

en temps reel avec un modele de description de I'ionosphSre, compare 
d des valeurs vraies ; 

la figure 6 est un graphique illustrant des exemples d'erreurs dans 
la determination du modele ionospherique en temps reel de la double 
10 difference de "STEC" compare a deux seuils ionospheriques 

les figures 7A a 7D sont des graphiques illustrant des exemples 
d'erreurs (en metres) estimes dans differents composantes de 
navigation ; 

les figures 8A £ 8D sont des graphiques illustrant des exemples 
15 d'erreurs (en metres) reelles dans les memes composantes de 

navigation ; 

la figure 9 est un graphique illustrant un exemple de pourcentage 
de reussite sur la determination instantanee d'une double difference 
ionosph6rique suffisante pour resoudre toutes les ambiguTtes (avec un 
20 "TECU" meilleur que 0,26) pour des corrections de temps de latence 

dans une gamme de 1 a 30 secondes ; et 

la figure 10 est un graphique illustrant un exemple de 
pourcentage de reussite sur la determination instantanee d'une double 
difference ionospherique pour des pour des corrections de temps de 
25 latence allant jusqu'a 900 secondes. 

Dans ce qui suit, sans en limiter en quoi que ce soit la port&e, on se 
placera ci-apres dans le cadre de ['application preferee de I'invention, sauf 
mention contraire, c'est-a-dire dans le cas d'un proced6 de navigation par 
satellite tr6s precis, d'ondes radioelectriques d trois porteuses de frequences 
30 distinctes emises par ces satellites et des corrections en temps reel issues d'un 
module en temps reel de description ionospherique d'une region sondee par le 
rayonnement radiofrequence de satellites. 
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On va maintenant d6tailler les principals etapes du procede selon 
I'invention, en faisant appel pour les illustrer a des valeurs numeriques 
typiques. II comprend trois 6tapes de base, qui sont semblables 3 celles du 
proc6d6 "TCAR" precite et vont £tre explicitees ci-dessous : 

5 Etape 1 : 

Pour rSsoudre I'ambiguitS de "chemin" dit "extra -long" (ou "extra-wide", 
selon la terminologie anglo-saxonne : typiquement avec une longueur d'onde 
de 7,5 metres sur un ensemble de donnSes d'un exemple d'une experience 
donn6e), il est ajoute une combinaison de codes dits de "pseudo-longueurs". 

10 Parmi les causes d'erreurs, on peut citer le ph6nomene connu dit de 

"trajets multiples". Bien que les trajets multiples de pseudo-longueurs puissent 
diminuer le taux de rSussite, cette erreur est en general de faible amplitude 
comparee a la grande longueur d'onde du "chemin extra-long", et il est 
habituellement possible de surmonter ce probleme. 

15 En effet, la premiere 6tape typique de I'approche "TCAR" (voir par 

exemple I'article precite de VOLLATH et al„ de 1998) est d'estimer la double 
difference d'ambigufte (entre des paires de recepteurs et de satellites), que Ton 
appellera ci-apres VAAT W , de la double difference de phases de porteuse de 
chemin "extra-long", que Ton appellera ci-apres VAL W (Lew 6tant la longueur 

20 d'onde de la porteuse de chemin "extra-long") en utilisant une combinaison de 
"pseudo-distances" (ou codes) Pew qui partagent les memes valeur et signe de 
la dependance ionospherique que les combinaisons de phase. Ceci est 
possible compte tenu de la grande longueur d'onde resultant de la 
combinaison de "chemin extra-long". Les valeurs de Lew et Pew peuvent etre 

25 determinees par les relations suivantes : 

P a> =^^=P+«^+M ew+ E a , + ... (2); 

relations dans lesquels L x est I'observation de phase de porteuse (en 
30 unites de longueur) a la frequence f x et a la longueur d'onde X*. Les valeurs de 
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X sont portees dans le TABLEAU I annexe en fin de !a presents description, 
etant entendu que les chiffres 1 § 3 sont assoctes aux trois frequences de 
Texemple decrit. Dans Le TABLEAU I sont rassembl6s un certain nombre de 
donnees associees : frequences, longueurs d'onde, erreurs de d6pendance 

5 ionospherique, de trajets multiples maximum, et d'observations diverses, utiles 
d la bonne comprehension du precede selon ('invention. Ces donnees sont 
associees aux trois bandes de frequences, L 1 a L 3 , respectivement, et aux 
phases de porteuses des chemins "extra-longs" et longs" (Lew et L w ) t 
respectivement. P x est I'observation correspondante de "pseudo-distance". Les 

10 erreurs de trajets multiples et d'observation pour des phases de porteuses et 
des "pseudo-distances" sont notees m Xi M Xi e x , E Xt respectivement (les valeurs 
maximaies de trajets multiples et les erreurs de mesure typiques sont 
egalement reportees dans le tableau I). L'ambiguite de phase de porteuse 
indiff6renciee Ax*>x, qui contient des retards instrumentaux et, apres double 

15 differentiation, devient une valeur entiere de multiplication de la longueur 
d'onde, X X V6N X . Le retard ionospherique aj est proportionel au "Contenu 
Total d'Electrons Obliques" ou "STEC", I'integration de la densite d'electrons 
libres le long du faisceau rayonne par un satellite etant habituellement mesure 
en unites TECU, comme il a et§ rappele. Dans les relations (1) et (2), p 

20 repr6sente un terme qui ne depend pas de la frequence (distance, erreurs 
tfhorloge, refraction troposphSrique, etc.). 

A partir de la relation (2), dans laquelle des termes additionnels mineurs, 
tels que le renroulement de phase de porteuse ne sont pas explicitement 
representee il est possible d'estimer la double difference d'ambiguTte de 

25 chemin "extra-long", que Ton appellera ci-apres VUf m , sur une seule "epoque" 
d'observation, en soustrayant le code correspondent, comme exprime dans la 
relation ci-dessous : 

VA# W ~VA(L W -Pj= VAtf w -^(VAM W + VAE W )+... (3) ; 



30 Quand il n'y a pas de trajets multiples importants affectant les recepteurs 

(usuellement inferieurs a 3,7 metres), la limite d'erreur d'une telle evaluation 
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est inferieure & 0,5 cycles (comme le montrent les valeurs foumies dans le 
TABLEAU I), ce qui permet de fixer instantan6ment Pambigurte a une valeur de 
nombre entier exacte. 
Etape 2 : 

5 L'ambiguTte de combinaison de "chemin long" est estimee a partir de la 

phase de porteuse non ambigue de "chemin extra-long", obtenue lors de la 
premiere etape. La difference entre celles-ci est constitute la plupart du temps 
par I'ambigurte de "chemin extra-long", et la refraction ionospherique 
differentielle (environ 0,06 cycles/TECU avec les frequences de travail du 

10 present exemple). Les termes non-dispersifs sont supprimes. Les problemes 
principaux sont ici I'erreur de mesure et les trajets multiples des signaux de 
phase de porteuse. Bien que les valeurs typiques de la refraction 
ionospherique differentielle aux latitudes moyennes et de lignes de base en 
dessous de 100 kilometres, soient justes de quelques TECU, I'utilisation de la 

15 correction ionospherique peut augmenter de maniere significative le 
pourcentage de reussite sur des distances plus importantes, et dans des 
scenarios ionospheriques plus difficiles. 

Les details principaux de cette etape vont maintenant etre explicites. 
Une fois que les premieres ambiguTtes de "grande longueur d'onde" 

20 VAAT W sont resolues en utilisant des "pseudo-distances" correspondantes, une 
deuxieme ambiguite telle que la combinaison de "chemin long" L w exprimee par 
la relation suivante : 

K = KdK -« 88 fA rJf 2 = p + KK +<*J + »*» + e w + ... (4) ; 

(relation dans laquelle fa et sont les phases de frequences U et L 2 .) peut 

25 etre estimee a partir de VAAT W et de la difference correspondante des phases 

de porteuse (voir les valeurs constantes correspondantes dans le TABLEAU I), 

conformement a la relation : 

VAN VAL w -VAL CT +A w VAAr w „ 
K 

VAN„ - ± VA(e w + m w - m w ) + ±{a w - a m )VA7 (5) ; 
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En effet, en presence de trajets multiples moderes, le terme d'erreur 
correspondent de la relation (5), qui inclut egalement I'erreur d'observation, est 
habituellement inferieure a 0,3 metres, c'est-a-dire inferieur a 0,4 cycles. Le 
terme d'erreur subsistant de la resolution d'ambiguYte correspond a la refraction 
5 ionospherique a une valeur de 0,0580 cycles/TECU pour les frequences de 
I'exemple (voir TABLEAU I). Ce terme peut degrader le taux de reussite, mais 
ce n'est pas un terme critique aux latitudes moyennes et aux distances 
inferieures a quelques centaines de kilometres, pour lesquelles la double 
difference des valeurs de "STEC", que Ton appellera ci~apr6s, VA/, sont en 

10 general inferieures a 1 0 TECU. Un exemple sera detaille ci-apres par reference 
au diagramme de la figure 5. 

C'est dans un tel contexte que ('utilisation d'un modele ionospherique 
en temps reel augmente les conditions de r6ussite de la resolution d'ambiguYte 
de "chemin extra-long", aux basses latitudes et pour des scenarios de lignes 

15 plus longues. 

foape 3 : 

L'ambiguYte de phase de U est derivee de la difference entre U et le 
"chemin extra-long" non ambigu obtenu precedemment. Lors de cette etape, le 
probleme principal est la r6fraction diff6rentielle ionospherique correspondante 
20 (environ 1,9 cycles/TECU), qui peut produire des erreurs de plusieurs cycles 
aux latitudes moyennes. 

Lors de la troisieme etape, une approche semblable a la deuxfeme 
etape est mise en oeuvre, mais en utilisant les differences de phase de 
porteuse entre de longueurs d'onde courte et moyenne, au lieu de longueurs 
25 d'onde moyenne et longue, comme il ressort de la relation suivante (6) : 

relation dans laquelle ct\ et ow (de fagon generate a x ) sont des coefficients 
ionospheriques definis dans le TABLEAU I, dans Pexemple decrit 

Lors de la troisieme etape, la combinaison de I'erreur de mesure de 
30 phase de porteuse et d'un trajet multiple modere introduit encore une erreur 
typique inferieure a 0,2 cycles (voir TABLEAU I). Cependant, le probleme 
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critique ici est la refraction ionospherique qui peut introduire des erreurs 
superieures a 0,5 cycles (-1,9475 cycles/TECU) egalement pour des lignes de 
base courtes. 

Cette limitation fondamentate peut £tre surmontee en calculant des 
5 corrections ionospheriques en temps r6el meilleures que 0,26 TECU (c'est-£- 
dire (0,5 cyclesx1TECU)/1,9475) pour assurer la bonne Evaluation d'ambigurte 
de nombre entier. Ces corrections sont introduites lors de P6tape 4. qui va 
maintenant §tre d6crite. 
El-tape 4 : 

10 Pour surmonter le probleme apparu a la troisieme 6tape (c'est-£-dire 

pour determiner les ambiguTtes de longueur d'onde les plus courtes), on 
determine un module en temps reel pour estimer la refraction differentielle 
ionospherique. Ce module est calcule a partir de donnees de phases de 
porteuses, a double frequence, £ des emplacements fixes d'un reseau, de 

15 stations de reference. II est assume qu'ainsi, il est realise une description de la 
region ionospherique traversee par le rayonnement d'ondes radioelectrique 
emises par des satellites. Les donnees de ce module sont combines de fagon 
connue avec des donn6es issues d'un calcul geodesique simultan6, ce dernier 
pouvant etre effectu6 avantageusement par Tune de stations terrestres fixes du 

20 reseau, dite maTtre. Lavantage principal de cette technique est qu'elle permet 
d'estimer la refraction differentielle avec une erreur typiquement inferieure a 
0,25 TECU, ce £ des distances de centaines de kilometres du site de reference 
le plus proche et dans des conditions ionospheriques variees. Cette precision 
implique des erreurs typiques en dessous de 0,5 cycles dans la fixation 

25 instantanee de PambiguTte U pour des distances moyennes et grandes. 

La distribution ionospherique d'&lectrons libres peut etre determinee 
de fagon approchee par une grille d'unites de resolution de volumes, ou 
"voxels", dans laquelle la densite d'electrons est assumee constante a un 
instant donne dans un systeme "Inertiel a Centrage Terrestre" ("Earth Centred 

30 Inertial" ou "ECI" selon Pabr6viation anglo-saxonne couramment utilisee). 

Une disposition typique de ce type est d6crite sur la figure 1A. Cette 
figure 1 A illustre schematiquement une tranche meridienne de voxels Vbx^ (/,/, 
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/c)etant des indices de coordonnees, pour la longitude, la latitude et i'altitude 
solaire, respectivement), dans lesquels la distribution ionospherique de 
densite d'electrons est decompose* selon la relation (7), de maniere a 
determiner les donnees d'un modele ionospheriques en temps reel. 

5 Sur cette figure 1A, on a represents le globe terrestre, GT, en coupe 

partielle, entoure de la couche ionospherique Con, couche que I'on a subdivise 
arbitrairement en deux couches intermediates, C« et C/ 2 , respectivement. 
L'altitude basse de la couche Q, est 60 Km, dans I'exemple decrit, et son 
altitude haute 740 Km. L'altitude basse de la couche C R est 740 Km et son 

10 altitude haute 1420 km. Les surfaces angulaires des voxels Vox^ sont de 5x2 
degrSs chacune. 

Pour fixer les idees, on a represents des satellites, SAT, a SAT n (trois 
ont ete representes sur les figures)dans lesquels sont embarquees des 
emetteurs "GPS", GPS* a GPSa, respectivement. Ces satellites, SAT, a SAT„, 

15 sont en vue de tout ou partie de stations terrestres. par exemple les trois 
stations terrestres fixes de reference, ST, a ST M , representees sur la figure 
1A, comprenant des receptees "GPS", GPS* a GPSr„ et des moyens de 
calculs integres (non representes). 

Sur la figure 1B, on a represents, de facon schematique, I'architecture 

20 d'un systeme de navigation complet, conforme a I'invention, vue du cdte du 
nomade SUR. 

Ce systeme n§cessite un Squipement, des interfaces et des 
ressources physiques, similaires a ceux de certains systemes actuellement 
dSployes, tels que EGNOS" ou WAS". 

25 Le systeme pour la mise en oeuvre du procede selon I'invention 

dernier comprend un recepteur "GPS", SURgps, qui recoit les signaux emis par 
des satellites, par exemple les satellites, SAT, a SAT n , de la figure 1A. A partir 
de ces signaux sur trois bandes de frequences, on effectue les trois Stapes du 
procSde TCAR", comme indique precedemment, a I'aide de moyens de calculs 

30 integres (non reprSsentes). 

II recoit egalement les signaux de donnees du modele ionospherique 
en temps reel determine par la station de reference REF (la station la plus 
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proche parmi les stations du r6seau de stations terrestres fixes de reference, 
S7*i a ST M , de la figure 1A) a partir de signaux re?us des satellites, S47V 
GPS E1 a SAT n -GPS Eni A I'aide un recepteur "GPS", REF GPSt . La station REF 
emet les donnees du modele vers le nomade SUR par I'intermediaire d'un 
5 6metteur REF Ej . 

Des donn6es de configuration g§od6sique, qui peuvent etre 
determines dans une seule station, dite maitre, REF Ml sont 6mis par son 
emetteur REF ME et distribues au nomade SUR, pour etre combines avec les 
donnees du modules ionospherique. Cette station comprend a priori egalement 
10 un recepteur "GPS", REFMgps, comme les autres stations terrestres fixes de 
reference. Les donnees de configuration geodesique sont calculees 
simultan§ment aux donnees du modele geod6sique. 

De fapon bien connue, la determination ionospherique peut etre 
obtenue en resolvant en temps reel, a Taide d'un filtre de Kalman la densite 
15 moyenne d^lectrons N e de chaque cellule illuminSe /, j, et k f (/, j et k etant des 
indices de coordonnees en trois dimensions, comme precedemment defini) 
traitee comme un processus de propagation al§atoire, et avec un bruit de 
processus typique de 10 9 a 10 10 electrons/m3/V^e«re . A titre d f exemple non 
limitatif, un un traitement d I'ide d'un filtre de Kalman est d6crit dans le livre de 
20 G.J. BIERMAN, intitul6 : "Factorization Methods for Discrete Sequential 
Estimation", paru dans "Mathematics in Science and Engineering", vol 128, 
"Academic Press", NEW YORK, 1977. 

Selon une caractfiristique importante du proced6 de Pinvention, 
seules les donnees de phase de porteuse sont utilises. De ce fait, le bruit de 
25 code de "pseudo-distance" et les trajets multiples sont evites. Les polarisations 
de phase de porteuse B f (constantes dans chaque arc continu de donnees de 
phase de porteuse pour chaque paire satellite-recepteur) sont estimees 
simultanement, en tant que variables aleatoires (qui deviennent des processus 
aleatoires de bruit blanc quand un glissement de cycle se produit). Dans le 
30 filtre, les polarisations sont decorrelees en temps r6el des valeurs de densite 
d^lectrons, dans la mesure ou la g6om6trie de satellite change et les 
variations des deux sortes d'inconnues sont devenues plus petites, comme le 



WO 2004/057364 



PCT/FR2003/050176 



18 

montre la relation (7) ci-apres qui repr6sente le modele pour des informations 
ionospheriques donnees, entre une paire satellite ["GNSS" - recepteur], U et 
L 2 6tant es phases de porteuse exprim6es en unites de longueur, Lf L r L 2 et 
N e la densite d'electrons). L t est donn§ par la relation suivante : 

L, = STEC + B t = C =EEE W Km + *' < 7) : 

JREC l j k 

5 relation dans laquelle les valeurs "REC" et "SAT 1 sont relatives a un recepteur 
GPS et un satellite, respectivement, d'une paire "recepteur-satellite", dl une 
distance elementaire sur la ligne de vis6e reliant cette paire et As UJ[ une 
surface elementaire du "STEC" de coordonnee /, j et k 6tant des indices de 
coordonnees en trois dimensions, comme precedemment defini. 
10 Cette approche convient, en particulier, pour detecter les 

caracteristiques locales de la distribution de density d'electrons, et ['utilisation 
de deux couches avec des donnees au sol d'un systeme de navigation mondial 
par satellite ("GNSS" selon I'abreviation anglo-saxonne) au lieu d'une seule 
couche, selon les precedes habituels, reduit de maniere significative le risque 
15 d'une mauvaise moderation du contenu d'electrons. 

Dans le cas des rSseaux "WADGNSS" (pour "Wide Area Differential 
GNSS" ou "GNSS" differentiels sur une large zone), b partir du contenu de ces 
corrections en temps reel du contenu total d'electrons obliques, ou "STEC", 
obtenues par la relation (7), il est possible : 
20 - preincrement, de former la double difference station-satellite, VASTEC, 

avec une erreur inferieure a 1 TECU et d'obtenir une deuxfeme ambiguTte 
(le "chemin long") dans les stations de reference ; et 
- deuxiemement, d'interpoler, dans le recepteur du nomade, la valeur U 
non ambigue, c'est-a-dire une valeur tres precise de V/SSTEC & un niveau 
25 de quelques centimes de TECU, valeur obtenue apres la resolution des 

ambiguTtes. 

Si la valeur interpolee est meilleure que 0,26 TECU, alors le recepteur 
du nomade peut rSsoudre les deux ambiguTtes en temps reel. 

Ce precede, en soi, est semblable au precede dit "WARTK" paru dans 
30 I'article de Hern^ndez-Pajares et al., (2000) pr6cite. Les resultats obtenus 



WO 2004/057364 



PCT/FR2003/050176 



19 

jusqu'ici avec la technique "WARTK", dans differentes experiences, sont 
recapitul6s dans un article egalement de Hern^ndez-Pajares et al., intitule : 
"Tomographic modeling of GNSS ionospheric corrections: Assessment and 
real-time applications", paru dans "ION GPS'2001", septembre 2001. On se 
5 reportera avec profit a ces deux articles pour une description plus detailtee de 
ce procSde. 

Pour I'utilisation cinematique de ce procede, une des limites les plus 
importantes est I'existence dlrregularites ionospheriques locales, telles que 
des perturbations Ionospheriques de propagation ou "TID" (abreviation anglo- 
10 saxonne de "Traveling Ionospheric Disturbances"), qui peuvent produire de 
mauvais resultats en utilisant une interpolation lineaire des corrections 
ionospheriques entre les stations de reference du reseau. Les performances 
peuvent etre amelior6es en incorporant des donnees duales de frequence du 
nomade. 

15 D'autre part, dans le cas du calcul dans la station de reference, les 

longues distances et les gradients forts de contenu d'electrons peuvent aussi 
limiter les performances de la technique. Pour cette raison, selon une 
caracteristique du procede de Pinvention une procedure supplemental a ete 
d6veloppee pour inclure le cas des stations permanentes et/ou de reference 

20 separees par des milliers de kilometres, en utilisant un code dit de "chemins 
longs" lisse, ce de fagon a ameliorer le modele ionospherique et par 
consequent constituer une aide ^ la determination en temps reel d'ambiguite. 
Une telle procedure est decrite, par exemple, dans Particle de Hernandez- 
Pajares et al., intitule : "Improving the real-time ionospheric determination from 

25 GPS sites at Very Long Distances over the Equator*, paru dans "Journal of 
Geophysical Research", 2002. 

Cependant les methodes mises en oeuvre dans Part connu pour 
resoudre, en temps reel, Pambiguite dans le recepteur d'un nomade, faisant 
appel a un systeme de frequences duales du "GPS" exige une grande quantite 

30 de calculs geodesiques, afin de garantir la determination de position 
relativement bonne, au niveau de la vingtaine de centimetres, pour resoudre 
des combinaisons de polarisations flottantes libres ionospheriques, et alors de 
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pouvoir (en mettant en ceuvre le precede "WARTK" precite) resoudre les 
ambiguites en temps reel, ce apres un temps de convergence de I'ordre de 
quelques minutes. 

Uinvention permet de pallier les deficiences de cet art connu, 
notamment presente par le procede "WARTK". Elle permet Amelioration des 
mSthodes instantanees de resolution d'ambiguYtes de phase de trois-porteuse 
aux distances moyennes et longues entre stations de reference (allant de 
dizaines a des centaines de kilometres), et avec un minimum de calculs 
geodesique. ^amelioration principale reside dans la troisieme etape du 
procedS "TCAR" precite, en utilisant des corrections ionospheriques en temps 
reel fournies par un module ionospherique calcule sans interruption dans les 
stations de reference. 

Une amelioration additionnelle de la premiere etape precitee est 
egalement implements dans le proc6d§ de I'invention. Cette amelioration 
consiste en ('utilisation de trois codes de pseudo-distances distincts, appeles 
ci-apres Pi, P 2 , P3, au lieu d'un element de code unique, ce qui permet de 
diminuer Teffet des trajets multiples dans la resolution d'ambiguTte de "chemin 
extra-large" 

En outre, dans un mode de realisation prefere, un essai d'integrite est 
realist en utilisant les codes de pseudo-distances de "chemin large" et de 
bande de frequence U pour d&ecter des sauts dans revaluation d'ambiguYte 
associee & une erreur dans la resolution d'ambiguYte longueur d'onde la plus 
longue. En effet, lors des deuxieme et troisieme etapes precitees, d partir des 
relations (5) et (6), on constate ais§ment q'un cycle dans Terreur d^mbiguTte 
de "chemin extra-large" produit environ huit cycles d'erreurs de "chemin large", 
et un cycle de "chemin large" se traduit par environ quatre cycles tferreur L u Li 
etant la porteuse de plus petite longueur d*onde dans I'exemple decrit. 
Souvent, ces sauts devaluation d'ambiguYte sont suffisamment importants pour 
etre detectes etfiltres en utilisant les codes correspondants. 

Cette approche, selon Tinvention, permet de venir a bout directement 
des limitations principales rendant impossible une navigation d une erreur 
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typiquement inferieure 3 quelques centimetres sur de longues distances (plus 
de 100 kilometres). 

Le TABLEAU II dispose en fin de la pr6sente description resume les 
caracteristiques comparees des principaux proced6s de resolution d'ambiguite 
5 en temps reel de Tart connu ("TCAR", "ITCAR", WARTK") et conforme a 
Pinvention. 

Pour mieux illustrer les caracteristiques principales du precede selon 
Invention, on va maintenant decrire un exemple d'experimentation realise sur 
la base de plusieurs jeux de donnees fournis par des genSrateurs de signaux 

10 modifies embarques dans des satellites d'un systeme existant du type "GNSS" 
precite, avec simulation d'aSronefs, d'utilisateurs nomades de "surface" 
(vehicules terrestre par exemple) ou de sites fixes, Des conditions dynamiques 
et ionospheriques variees ont et6 considerees. 

Pour fixer les idSes, a titre d'exemple significatif, la resolution 

15 d'ambiguite en temps r6el du recepteur d'un vehicule de surface, appele ci- 
apres nomade SUR, a ete etudi6e en detail, ce par rapport a une station fixe de 
reference, appetee ci-apres REF, 6loignee d'environ 129 km du nomade SUR, 
ce dans differents scenarios. Dans Texemple decrit, deux porteuses de "GPS" 
(1575,42 et 1227,60 MHz respectivement), et la porteuse du canal 24 du 

20 systeme "GLONASS", de frequence 1615,50 MHz, constituaient les trois 
frequences adoptees par un simulateur du systeme "GNSS", pour quatre 
satellites en vue pendant 20 minutes a 1Hz, en utilisant les douze canaux 
disponibles d'un recepteur de validation de type connu sous la reference 
"AGGA" (pour "Advanced GPS/GLONASS ASIC"). "AGGA" est un circuit 

25 integre numerique permettant un traitement de signal numerique a haute 
Vitesse pour des applications spatiales, telles que les occultations de radio- 
transmissions d signaux de systemes "DNSS" et de d6termination d'orbites 
terrestres basses. Ce circuit est realise a base d'un "ASIC" standard (ou "circuit 
integre pour applications spScifique", dans Texemple d6crit, un composant 

30 T7905E de la societe ATMEL). "GLONASS" ("Global Orbiting Navigation 
Satellite System") est un systeme de navigation russe par satellite comparable 
au systeme "GPS". 
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Deux ensembles de donnees de base ont ete considers dans 
I'exemple decrit, ci-apres references "P5-M0", ensemble associe a une 
puissance maximale de signal et pas de trajets multiples, ce qui constitue un 
cas ideal, et "P3-M1", ensemble associe a une puissance moyenne et des 

5 trajets multiples. 

En plus des recepteurs des stations de reference, simules dans les 
ensembles de donnees, trois stations additionnelles, existant dans le reseau 
"IGS" (pour "International GPS Service" ou "Service International GPS"), ont 
ete ajoutees a des distances superieures a 200 km, uniquement pour le calcul 
10 ionospherique, ce qui a permis d'emuler une situation plus realiste comprenant 
un reseau de sites fixes plus important. 

On va maintenant montrer les resultats d'une etude detaillee pour le 
cas plus difficile, c'est-a-dire la ligne de base de la paire "nomade SUR - 
station REF' la plus longue, soit 129 km environ. La station de reference REF 
15 est I'une des stations ST, a ST M de la figure 1 A. 

La figure 2 illustre cet exemple et represente la ligne de base et les 
positions de la station REF et du nomade SUR, par rapport a des axes de 
coordonnees, latitude ellipsoTdale et longitude, gradues en degres. 

La figure 3 est un agrandissement (zoom) de la figure 2 montrant la 
20 trajectoire du nomade SUR, toujours par rapport a des axes de coordonnees, 
latitude ellipsoTdale et longitude, gradues en degres. 

Le premier aspect des resultats est reference en tant que performance 
du filtre ionospherique en temps reel. Des exemples de retards estimes du 
modele ionospherique en temps reel, references Tm sur le graphique, 
25 compares aux retards ionospheriques reels, references W, sont represents 
sur la figure 4. Ces valeurs correspondent a des mesures d'une 
experimentation effectuee pour montrer les possibility du procede selon 
I'invention. 

Et, parametres plus importants pour une navigation precise, les 
30 valeurs et les differences (c'est-a-dire les erreurs) entre les doubles differences 
de "STEC" estimees, VA/, pourle nomade SUR, et leurs valeurs reelles, VA/ , 
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sont reportees sur les figures 5 et 6 (en ordonnee, exprimees en TECU ; 
Tabscisse representant des retards en secondes). 

De fagon plus precise, la figure 5 illustre la double difference de 
"STEC", notee n ddSTEC\ estimee en temps reel avec un modele 
ionospherique, conformSment d Tun des aspects importants de Invention, 
compart aux valeurs vraies pour le recepteur nomade SUR, ce rapporte au 
rScepteur le plus eioigne, soit REF (eloignement d'environ 129 km). 

La figure 6 illustre I'erreur dans la determination du modele 
ionospherique en temps r§el de la double difference de "STEC", 6galement 
notee "ddSTEC". Les valeurs reelles du modele ionospherique sont comparees 
a deux seuils ionospheriques, "Seu//i" (+ 0,26 TECU) et "Seuil 2 tt (- 0,26 TECU), 
permettant de determiner la troisieme ambiguite VAtf, . 

Plus de 92 % des determinations ionospheriques pr6sentent des 
erreurs au-dessous de la limite de seuil de ± 0,26 TECU. Elles sont done assez 
precises pour pouvoir r6soudre les trois ambigurtes en Tabsence de trajets 
multiples et pour ne pas tenir compte de I'erreur de mesure. La majeure partie 
des 8 % devaluations restantes presentant des erreurs superieures & 
0.26TECU provenaient d'un satellite qui etait observe a basse altitude et dans 
la direction du sud, direction pour laquelle les gradients ionospheriques sont 
les plus eleves. Ces r6sultats ont pu encore etre legerement ameliores en 
additionnant aux quatre stations de reference une cinquieme effectuant des 
calculs ionospheriques de corrections pour la direction sud. 

Une fois que les corrections ionospheriques en temps reel ont 6t£ 
calculees et transmises a partir du reseau de reference, les Stapes du procede 
selon Tinvention peuvent etre effectuees dans le nomade SUR, de la fagon 
precedemment decrite pour estimer et resoudre instantanement ("single 
epoch") les trois ambiguTtes par faisceau. 

Un resume des principaux resultats de la resolution d'ambigu'fte, pour 
les ensembles precites de donnees "P5-M0" et "P3-M1" respectivement, est 
presente sous la forme de deux TABLEAUX, III et IV, disposes en fin de la 
presente description. Les taux de reussite des trois 6tapes du proc6d6 "TCAR", 
sont indiques dans les trois cas suivants : 
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(a) sans corrections ionospheriques ; 

(b) avec des corrections ionosph6riques correspondantes d'un 
modele de corrections ionospheriques de Klobuchar, emises par le 
systeme "GPS " actuel ; et 

5 (c) avec des corrections de module ionosph6rique en temps r6el, 

conform6ment au procede de I'invention. 

De fagon plus precise, dans Pexemple d6crit en relation avec le 
TABLEAU III, le taux de reussite (en %) est relatif d un nombre de 3834 essais 
de resolution au vol pour les ambiguTtes de "chemin extra-long", de "chemin 

10 long" et de frequence de la bande Li (respectivement VAW W , VAAf w et VzW,) 
pour le r6cepteur nomade SUR, rapporte au site fixe plus eloign^ REF 
(eloignement d'environ 129 kilometres), et avec I'ensemble de donnees id6al 
P5-M0 precite (puissance maximale de reception et aucun trajet multiple). 

Le TABLEAU IV est similaire au TABLEAU III, mais pour I'ensemble de 

15 donnees P3-M1 precite (niveau de puissance de reception moyen et trajets 
multiples). Les chiffres entre parentheses indiquent le taux de reussite ramene 
au nombre total d'observations. 

On constate aisement, a I'examen des TABLEAUX III ET IV, que le 
proced6 de Invention combinant les trois Stapes du procede "TCAR" avec une 

20 etape duplication d'un modele ionosph6rique precis en temps reel, permet 
d'augmenter considerablement le taux de reussite de determination d'ambiguite 
instantanee ("single epoch"), non seulement dans un scenario ideal (ensemble 
de donnees P5-M0 : de 0 d 92 %) mats egalement dans un scenario plus 
difficile (ensemble de donnees "P3-M1") avec des trajets multiples (35 % meme 

25 en utilisant une correction du type "modele de Klobuchar", alors que taux de 
reussite atteint 92% pour VAiV^ avec des corrections du modele 

ionospherique realisees selon le procede de Pinvention). 

L'examen du TABLEAU IV montre qu'un des problemes importants 
avec les ensembles de donnees "P3-M1", affectes par des trajets multiples, est 
30 le manque d'integrite dans les evaluations d'ambiguTtes selon le procede 
"TCAR", VAW W et VAtf w , avec environ 10 % et 4 % de determinations 
instantanees incorrectes, respectivement. 
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Pour augmenter Tintegrite, c'est-a-dire pour diminuer la possibility de 
prendre en consideration des ambiguites incorrectes, dans ce cas-ci en 
presence de phase de porteuse et pour coder par trajets multiples, il est 
possible d'utiliser les codes correspondants de pseudo-distances de "chemin 

5 long" et L u pour tenter de detecter des erreurs potentielles de determination 
d'ambiguTte de Pambiguite de longueur d'onde la plus longue precedente dans 
I'approche "TCAR" (des ambiguites de "chemin extra-long" et de "chemin long", 
respectivement), erreurs qui sont amplifies par un facteur de 9 et 4 longueurs 
d'onde dans le "chemin extra-long" et le "chemin long", respectivement, comme 

10 le montre les relations (5) et (6). 

Un resume des resultats correspondants est presente sous la forme de 
deux TABLEAUX, V et VI, disposes en fin de la presente description, en 
utilisant respectivement les pseudo-distances et les pseudo-distances lissees 
pour filtrer de telles erreurs importantes. 

15 Le TABLEAU V est similaire au TABLEAU II, mais le taux de reussite 

de resolution d'ambigu'ite est ici calcute apr6s avoir subi un test d'integrite, 
c'est-a-dire une comparison a I'ambiguYte d6riv6e de pseudo-distance pour 
filtrer les grands sauts dans 1'ambiguTte de "chemin long" et de U ces sauts 
etant dus a des erreurs pr6c6dentes d'ambiguTte de longueur d'onde plus 

20 importante. La dernfere colonne indique la disponibilite, c'est-a-dire le 
pourcentage des 3834 observations passant I'essai d'integrite precedent avec 
les pseudo-distances. Le signe entre parenthese "*" indique les pourcentages 
des ambiguites calculees apres passage du test d'integrite. Comme 
pr6c6demment, les chiffres entre parentheses indiquent le taux de reussite 

25 ramene au nombre total d'observations. 

Le TABLEAU VI est similaire au TABLEAU III, mais des pseudo- 
distances lissees sont utilisees a la place de pseudo-distances "pures" pour 
realiser le test d f integrit6. 

On constate aisement que procede de Tinvention permet une 

30 amelioration significative d'environ 20 % de I'integrite (de 79 % : TABLEAU IV a 
91 % : TABLEAU VI), ce avec une diminution relativement faible (16 %) de la 
disponibilite (celle-ci passe de 100 % a 84 %) par I'utilisation d'un test 
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d'integrite de code Iiss6. Si des codes instantanes ("single epoch") sont utilises 
en remplacement, I'integrite est 6ga!ement ameliorate (84 %) et, ce, avec une 
disponibilite plus elev§e (90 %). En utilisant un modele connu d'emission 
"GPS" de type Klobuchar, au lieu des corrections du module ionospherique en 

5 temps r6el selon le procede de I'invention, on constate une deterioration nette 
des resultats, accompagne d'un manque quasi-complet de disponibilite. 

Les ensembles de donnees disponibles pour I'experimentation qui 
vient d'etre detaillee, malgre le fait qu'ils etaient approprtes pour etablir une 
resolution instantan6e d'ambiguTte, etaient cependant tres limites pour 

10 determiner une position de fapon instantanSe ("single epoch"). En effet, pour 
('experimentation, seuls quatre satellites en vue etaient disponibles pendant 20 
minutes de donnees, de telle sorte que la solution de navigation instantanee 
pour le nomade SUR a du etre calculee en utilisant I'ensemble de donnees 
"P5-M0" sans trajets multiples, et avec un rapport signal sur bruit disponible 

15 tres eleve, afin de disposer rfun nombre minimum de quatre satellites pour la 
navigation. De plus, le retard tropospherique, pour revaluation duquel un 
cinquieme satellite au moins aurait 6t§ nScessaire, a ete omis. Enfin, le premier 
intervalle, entre 7200 a 7500 secondes approximativement, a ete evite dans I 
determination du positionnement en relation avec certaines configurations de 

20 coupure de phase de porteuse, potentiellement dus a des problemes de 
mesure dans le recepteur pendant cette periode. D'autre part, et du point de 
vue des corrections ionosph§riques et du taux de r6ussite correspondant de la 
resolution d'ambiguTte, les r6sultats sont pratiquement les memes que ceux 
obtenus pour Tensemble de donnees "P3-M0" (voir TABLEAU III). 

25 Ces limitations sont dues a I'utilisation d'ensembles de donnees 

foumis par des systemes existants, n'ayant subi qu'un minimum rfadaptations. 
Avec un ensemble de donnees plus realiste, par exemple fourni par six 
satellites ou plus, ces limitations importantes concernant la determination de 
positionnement disparaTtront. II sera en effet possible d'obtenir une meilleure 

30 g6orrtetrie (un facteur de geometric "DOP", ou facteur de "Dilution de 
Precision" de geometrie, plus faible) et la possibilite de detecter et de filtrer des 
satellites (typiquement aucun, ou parfois un seul) presentant des erreurs 
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d'ambiguTte. Ceci pourra etre realise en utilisant un filtre de navigation ou des 
algorithmes de surveillance autonomes de recepteur). 

Les r6sultats correspondants aux quatre satellites contenus dans 
I'ensemble de donnees P5-M0 sont indiques dans les figures 7A a 7D, 

5 montrant les r6siduels de pre-ajustement (Figure 7A) f les composants d'Est 
(Figure 7B), de Nord (7C) et de verticale (7D) de I'erreur de positionnement 
instantanee apres mise en ceuvre du precede selon I'invention, pour un 
nomade SUR eloigne de 129 km environ d'une station de reference REF. les 
ordonnees des graphiques sont graduees en metres et les abscisses (retards) 

10 en secondes. 

A titre de comparaison, les figures, 8A a 8D, represented des 
donnees correspondantes mais pour des ambiguTtes reelles. 

La caracteristique principale des residuels de pre-ajustement VAL C 
(Figure 7A) sont des erreurs d'environ 10 cm assoctees a une erreur de un 

15 cycle dans U (environ 8 % : voir TABLEAU III), qui produisent des erreurs de 
navigation ("single epoch") semblables a des sauts. Elles sont egalement 
representees dans les figures 7B a 7D. Elles sont amplifies par la dilution 
correspondante de facteur geom6trique de precision ("DOP"). Ces mauvaises 
Evaluations de TambiguTte U affectent typiquement une double difference sur 

20 chacune des trois qui sont disponibles par epoque (le nombre minimum pour 
determiner un positionnement), ce de telle mantere que cette erreur affecte 
trois fois le positionnement (environ 24 % des epoques, avec 76 % des erreurs 
en trois dimensions en dessous de 5 cm et 100 % en dessous de 21 
centimetres). Bien que la distribution des erreurs ne soit pas gaussienne, les 

25 valeurs efficaces ("RMS") resultantes en trois dimensions sont de 7, 3,5 et 2 
centimetres, pour les composants d'Est, de Nord et de verticale, 
respectivement. 

Pour resumer ces resultats, des erreurs de 1, 1 et 2 cm "RMS" sont 
obtenues pour les composantes d'Est, de Nord et de verticale, respectivement, 
30 quand I'ambiguTte est bien resolue (environ de 92 % des essais et 77 % des 
6poques) et de 3, 5 et 2 cm "RMS" pour les composantes d'est, de Nord et de 
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verticale, respectivement, en incluant les epoques avec une mauvaise 
resolution d'ambigui'te. 

En outre, les observations minorees de calculs de r6siduels, avant 
ajustement, de la double difference de phase de porteuseset les erreurs de 

5 navigation instantanees ("single epoch") sont representees sur les figures 8A d 
8D. On constate que I'amplification du bruit de phase de porteuse est nette, en 
particulier dans la composante verticale, atteignant environ 5 cm d'erreur. 
Cette trajectoire (avec une valeur "RMS" en trois dimensions de 3 cm, 1,1 cm 
et 2 cm pour les composantes de Nord, d'Est et de verticale, et des niveaux de 

10 95 % a environ 2, 2 et 4 cm, respectivement) represente tout a fait bien la 
solution qui pourrait etre obtenue avec un filtre en temps reel, au lieu de la 
solution en mode instantane ("single epoch"), apres la resolution correcte des 
ambiguTtes pendant les premieres epoques. 

Afin de caracteriser I'impact des temps de latence dans la correction 

15 ionospherique (dO par exemple a des problemes potentiels dans les 
communications), des retards de 1 a 30 secondes (30 Epoques) ont et6 
consideres pour les calculs ionospheriques de correction dans le reseau de 
stations fixes. Pour chacun de ces retards, on a calcule le pourcentage de 
reussite pour une "ddSTEC" (en TECU) pour obtenir une precision absolue de 

20 0,26 TECU, pour un nomade 6loign6 de 129 km. Ce pourcentage a <§t6 reports 
sur la figure 9 en fonction du temps de latence (en secondes). On constate 
aisement que Teffet des temps de latence est negiigeable jusqu'a 30 s, avec un 
taux de reussite atteint, pour tous les temps de latence, superieur a 90 %. 

Sur la figure 10, de plus grands temps de latence ont ete consideres, 

25 mais ceci a la fois pour les corrections de reseau de reference et de nomade, 
de telle sorte que les points d'entr6e ionospheriques restent etroits. On 
constate qu'apres 5 minutes, le taux de reussite diminue de 90 % a 85 %. 
AprSs 10 minutes, le taux de r6ussite diminue encore pour atteindre 75 %. Ces 
valeurs peuvent se degrader plus fortement dans des scenarios presentant une 

30 variation du retard ionosph6rique differential encore plus elevee. 
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On constate done, qu'a une latitude moyenne, le temps de latence ne 
pose pas un probieme important pour le proc6d6 selon I'invention qui peut 
supporter des temps de latence allant typiquement jusqu'a 5 minutes environ. 

A la lecture de ce qui precede, on constate aisement que I'invention 
5 atteint bien les buts qu'elle s'est fixes. 

Le procede selon I'invention permet, en appliquant des corrections 
ionospheriques en temps reel, £ I'aide d'un module tomographie de 
I'ionosphere, accompagn6es d'un minimum de calculs g6odesiques, une 
resolution instantanee d'ambiguite complete. Elle permet une navigation avec 
10 une precision de I'ordre du centimetre, £ des distances de plus de 100 
kilometres de stations de reference. 

L'invention trouvera une application pr6feree dans les prochains 
systemes de navigation par satellite, tel que "G ALU LEO" ou le "GPS" 
modernise, mettant en ceuvre remission de porteuses sur trois bandes 
15 frequences distinctes. 

Selon une caracteristique principale du procede selon I'invention, 
cette possibilite de navigation instantanee, a une precision meilleure que le 
decimetre, est due au fait que Ton combine un procede de resolution 
d'ambiguite de phase de trois porteuses, de type TCAR", peu complexe, avec 
20 une determination d'un modeie ionospherique en temps reel en vue de pouvoir 
appliquer des corrections ionospheriques. 

Les exemples d'experimentation rappeies dans la presente 
description, montrent montre une amelioration nette des resultats obtenus par 
des precedes de Part connu, dans des conditions equivalentes : 
25 - Le taux instantane de reussite C'single epoch") d'environ 60 % ou 

moins est am6liore et passe a environ a 90 %, ce e des distances de plus 
de 100 km d'une station terrestre fixe de reference la plus proche, et dans 
des conditions difficiles pour realiser un modeie ionospherique (& midi au 
maximum solaire) ; 

30 - Ce rendement 6leve est maintenu avec des temps de latence d'environ 

5 minutes dans les corrections ionospheriques ; 
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- La navigation correspondante, malgre le nombre limits de quatre 
satellites disponibles, pour les conditions d'experimentation retenues a 
titre d'exemple, dans les ensembles de donnees, peut etre fait 
instantan^ment, en mode "single epoch", et permet des valeurs "RMS" en 

5 trois dimensions, de 3, 5 et 2 centimetres pour les composantes d'Est, de 

Nord et de verticale, respectivement ; 

- Dans un mode de realisation pr6fere, un essai en temps r6el d'integrit§ 
d'ambiguTte, dans la mesure ou une approche de reduction simple de 
trajets multiple pseudo-distances a 6t6 incorpor§e au precede a permis 

10 d'en am&iorer les performances ; et 

- Ces bonnes performances ont 6te 6galement prouvSes dans des 
scenarios difficiles : 

. Faible latitude (35 degres au lieu de 48 degres), avec des valeurs 
ionospheriques plus elevees. 
15 . Conditions extremes en ce qui concerne le tropique, au-dessous de 

TAnomalie £quatoriale Nordique" de I'lonosphere, ou se produisent 
typiquement les gradients les plus eleves. Dans ce dernier cas, il est 
n§cessaire que Tutilisateur produise sa propre alimentation de filtre 
ionospherique et la combine avec les corrections du r6seau de 
20 reference et ses propres observations ionospheriques. 

. Dynamique elevee, rendant possible la fourniture de corrections 
ionospheriques precises utilisable par un aeronef volant d'avion a 
environ 140 km de la station de reference le plus proche. 
II doit etre clair cependant que invention n'est pas limitee aux seuls 
25 exemples de realisations explicitement decrits, notamment en relation avec les 
figures 1 a 8. 

Enfin, les exemples numeriques rfont ete fournis que pour mieux fixer les 
idees et ne sauraient constituer une quelconque limitation de la portee de 
invention, lis precedent d f un simple choix technologique a la portee de 
30 PHomme de Metier. 
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TABLEAU I 





X=1 


X=2 


X=3 


X=ew 


X=w 


Frequence 

MHZ) 


1575.42 


1227.6 


1615.5 


40.08 


347.82 


Longueur d'onde 

{X x , m) 


0.1903 


0.2442 


0.1856 


7.4799 


0.8619 


Erreur de mesure de 
phase 

(£^,171) 


r\ ftno 


«U.UUZ 






~U.U I 


Erreur de mesure de 
code 

(E^m) 


=3 


=3 


=3 


«2.1 


=2.1 


Phase de trajets 
multiples maximale 

(m xt m) 


«0.05 
(50.01) 


«0.06 
(50.01) 


«0.05 
(50.01) 


«2 
(50.5) 


«0.2 
(50.05) 


Code de trajets 
multiples maximum 

{M Xi m) 


«450 
(510) 


«450 
(510) 


«450 
(510) 


«450 
(510) 


«450 
(510) 


Coefficient 
ionospherique 

(a x ,mfTECU) 


-0.1623 


-0.2673 


-0.1543 


0.2083 


0.2031 



TABLEAU II 





Avantages 


Inconvenients 


"TCAR" 


Faible charge de calculs 


Erreur ionospherique limitant 
serieusement la resolution de 
troisieme ambigurte 


"ITCAR" 


Resultats ameliores en integrant 
TCAR dans un filtre de navigation 


Le retard ionospherique 
continue de limiter la resolution 
de troisieme ambiguTte 


"WARTK" 


Une modelisation en temps r6el 
precise permet une navigation 
precise a une centaine de 
kilometres de la station de 
reference la plus proche 


Charge elevee de calculs : un 
besoin de calcul d'une 
premiere solution independante 
de I'ionosphere pour le nomade 
et un temps de convergence 
significatif sont necessaires 


Proc6de 
conforme 
a 

I'invention 


Faible charge de calcul et 
modelisation precise en temps 
reel de Pionosphere, permettant 
une navigation du type "single 
epoch" precise a des 
gloignements de plus de 100 Km 
d'une station de reference. 
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TABLEAU III 



P5-M0 / SUR2-REF5 
(=129km) 


Taux de 
reussite 


Taux de 
reussite 


Taux de 
reussite 

VAW, 


Sans corrections 
ionospheriques 


100% 


100% 


0% 


Corrections de Klobuchar 


100% 


100% 


33% 


Corrections 
ionosoheriques 
en temps reel 


100% 


100% 


92% 



TABLEAU IV 



P5-M0 / SUR2-REF5 

(<=129km) 


Taux de 
reussite 


Taux de 
reussite 


Taux de 
reussite 

VAAT, 


Sans corrections 
Ionospheriques 


90% 


95% (86)% 


3% (2) 


Corrections de Klobuchar 


90% 


95% (87)% 


35% (31) 


Corrections 
ionospheriques 
en temps reel 


90% 


95% (86)% 


92% (79) 



TABLEAU V 



P3-M1 / SUR2-REF5 
(«129km) 


Taux de 
reussite 

VAN * 


Taux de 
reussite 

VAN W * 


Taux de 
reussite 

VAAT,* 


Disponibilite 


Sans corrections 
ionospheriques 


99% 


95% (94) 


0(2)% 


38% 


Corrections de 
Klobuchar 


99% 


97% (96) 


33% (32) 


38% 


Corrections 
ionosoheriques 
en temps reel 


96% 


96% (92) 


91% (84) 


90% 



5 TABLEAU VI 



P3-M1 / SUR2-REF5 

(=129km) 


Taux de 
reussite 

VAN m * 


Taux de 
reussite 

VAN„* 


Taux de 
reussite 

VAN, 


Disponibilite 


Sans corrections 
ionospheriques 


100% 


100% (100) 


0% (0) 


0.002% 


Corrections de 
Klobuchar 


100% 


100% (100) 


37% (37) 


0.002% 


Corrections 
ionosoheriques 
en temps reel 


100% 


100% (100) 


91% (91) 


84% 
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REVEND1CATI0NS 

1, Precede de navigation en temps r6el & I'aide de signaux radioelectriques 
a trois porteuses de premiere, deuxfeme et troisieme frequences distinctes, 
de valeurs croissantes en allant de ladite premiere a ladite troisieme 
frequence, pour determiner la position d'un utilisateur, dit nomade, lesdits 
signaux radioelectriques etant emis par un nombre d6terrnin6 d'§metteurs 
embarques dans des satellites toumant autour du globe terrestre et en vue 
dudit nomade, lesdits signaux etant re?us par un r6cepteur associe au dit 
nomade et un recepteur associe a au moins une station terrestre parmi une 
pluralite de stations terrestres fixes, dites de reference, lesdits signaux 
radioelectriques traversant une couche de I'atmosphere dite ionosph6rique 
entourant ledit globe terrestre et y subissant des perturbations engendrant 
des ambiguites de phases desdites porteuses, caracterise en ce qu'il 
comprend au moins les etapes suivantes : 

une premiere etape consistant en la determination, dans ledit 
nomade (SI//?), de FambiguTte de difference de phases entre 
lesdites troisieme et deuxieme porteuses, dite de "chemin extra- 
long", & partir d'une combinaison de pseudo-distances utilisant une 
valeur de code unique ; 

une deuxieme 6tape consistant en I'estimation, dans ledit nomade 
(SUR), de Tambiguite de difference de phases entre lesdites 
premiere et deuxieme porteuses, dite de "chemin long", a partir de 
ladite ambiguite de "chemin extra-long" d§terminee lors de ladite 
premiere etape ; 

une troisieme 6tape consistant, dans ledit nomade (SUR), en la 
resolution de Tambiguite d'une desdites frequence, a partir de ladite 
ambiguite de "chemin long" estim6e lors de ladite deuxieme etape ; 
et 

une etape supplemental Tapplication de corrections 
ionospheriques en temps r£el pendant ladite troisieme §tape, 



WO 2004/057364 



PCT/FR2003/050176 



34 

lesdltes corrections ionospheriques etant basees sur un modele 
ionospherique en temps r6el de ladite couche ionospherique (C/ow) 
continument mis £ jour. 

2. Procede selon la revendication 1, caracterise en ce que lors de ladite 
5 troisieme etape, la dite resolution d'ambigurte est effectuee sur ladite 

premiere frequence de porteuse. 

3. Procede selon Tune quelconque des revendications 1 ou 2, caracterise en 
ce que ledit modele est un modele ionospherique de description de ladite 
couche ionospherique, determine par au moins une desdites stations 

10 terrestres de reference (REF) recevant des signaux emis par un nombre 
predetermine desdits satellites {SATrGPS E i a SAT^-GPSe*) toumant autour 
du globe terrestre (G7) et en vue de celle-ci, lesdits signaux comprenant au 
moins deux porteuses de frequences distinctes, en ce que ledit mod§Ie est 
determine a partir de donnees de phases desdits signaux emis, et en ce qu'il 

15 comprend une etape de transmission de donnees correspondent au dit 
modele ionospherique. 

4. Precede selon la revendication 3, caracterise en ce que ladite 
determination de modele ionospherique est obtenue a partir de I'estimation 
de la distribution d'eiectrons libre dans ladite couche ionosphere (C /0 w), en 

20 ce que cette estimation etant realisee de fagon approchee en decomposant 
la couche ionospherique (C /0 w) en une grille d'unites de resolution de 
volumes (Vox^), dits "voxels", illuminees par le rayonnement radioeiectrique 
desdits signaux se propageant dans la dite couche ionospherique (C /0 w), 
dans lesquelles la distribution ionospherique de densite ^electrons est 

25 assumee constante S un instant donne, et en ce que ladite determination est 
obtenue par resolution en temps r6el de la densite moyenne d'electrons 
dans chacune desdites unites de volume {Voxq k ) illuminees par ledit 
rayonnement radioeiectrique en faisant usage d'un filtre dit de Kalman. 

5. Procede selon ia revendication 4, caracterise en ce qu'il comprend une 
30 etape supplemental consistant d combiner des donnees associ6es au dit 
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moddle ionospherique a des donnees geodesiques calculees 
simultanement, et en ce que lesdites donnees geodesiques sont calculees 
par une seule des dites stations terrestres fixes de reference (REFm-REFme), 
dite station maTtre, et distributes a ladite pluralite de stations terrestres fixes 
5 de reference (REF). 

6. Proc6de selon la revendication 1, caracterise en ce qu'il comprend une 
etape supplementaire consistant a utiliser trois codes de pseudos-distances, 
associes aux dites trois porteuses, lors de ladite premiere etape 
determination de I'ambiguite de difference de phases entre lesdites 

10 troisieme et deuxieme porteuses, dite de "chemin extra-long". 

7. Proc6de selon la revendication 2, caracterise en ce qu'il comprend une 
etape supplementaire consistant a effectuer un essai d'integrite en utilisant 
des codes de pseudo-distances de tt chemin large" et de ladite deuxieme 
frequence pour detecter des sauts associes a une erreur dans ladite 

15 resolution d'ambiguite de ladite premiere frequence. 

8. Systeme de navigation par satellite pour la mise en ceuvre du procede 
selon Tune quelconque des revendications prec6dentes, caracterise en ce 
qu'il comprend une pluralite de satellites (SAT r GPS E , a SAT n -GPS En ) 
tournant autour du globe terrestre (GT), chacun des satellites emettant 

20 lesdits signaux a trois porteuses de frequences distinctes, au moins un 
nomade (SUR) comprenant un rScepteur (SUR GPS ) desdits signaux et des 
moyens de calcul integres effectuant lesdites premiere a troisieme etapes et 
integrer lesdites corrections ionospherique issues d'un modele 
ionospherique de description d'une region de I'ionosphere traversee par le 

25 rayonnement radioelectrique desdits signaux emis par ladite pluralite de 
satellites (SAT^GPS m a SAT^GPS^), une pluralite de stations terrestres 
fixes dites de reference (REF) comprenant chacune un r6cepteur (REFgps) 
recevant lesdits signaux emis par lesdits satellites, des moyens de calculs 
integres pour la determination dudit module ionospherique de description de 

30 la couche ionospherique (C/ow), et un emetteur (REF E ) pour transmettre des 
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donn6es correspondant au dit modele ionospherique au dit recepteur 
(SURgps) de nomade (SUR) et en ce qiTau moins une desdites stations 
terrestres fixes de reference (REF M ), dite maTtre, comprend un recepteur 
(REFMgps) desdits signaux emis par ladite pluralite de satellites (SAT r 
GPS E <\ a SAT n -GPSai), des moyens de calcul de donn6es geodesiques et un 
6metteur (REFme) pour les distribuer a ladite plurality de stations terrestres 
fixes de reference (REF). 

9. Systeme selon la revendication 8, caracterise en ce que ledit nomade 
(SUR) est situ6 d une distance de plus de 100 Km de la station terrestre fixe 
de reference la plus proche (REF). 
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